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Few examples of magnetic systems displaying a transition to pure dipolar magnetic order are 
known to date. As was recently shown, within the newly discovered class of single molecule magnets 
quite attractive examples of dipolar magnetism may be found. The molecular cluster spins and 
thus their dipolar interaction energy can be quite high, leading to reasonably accessible ordering 
temperatures even for sizable intercluster distances. In favorable cases bonding between clusters in 
the molecular crystal is by Van der Waals forces only, and no exchange paths of importance can be 
distinguished. An important restriction, however, is the requirement of sufficiently low crystal field 
anisotropy for the cluster spin, in order to prevent the occurrence of superparamagnetic blocking at 
temperatures above the dipolar ordering transition. This condition can be met for molecular clusters 
of sufficiently high symmetry, as for the Mn6 molecular cluster compound studied here. The uniaxial 
anisotropy of the cluster spin S = 12 is as small as D/ks ~ 0.013 K, giving a total zero-field splitting 
of the S — 12 multiplet of 1.9 K. As a result, the electron-spin lattice relaxation time remains fast 
(~ 10 -4 s) down to T c and no blocking occurs. Magnetic specific heat and susceptibility experiments 
show a transition to ferromagnetic dipolar order at T c = 0.16 K. Classical Monte Carlo calculations, 
performed for Ising S = 12 dipoles on a lattice do predict ferromagnetic ordering and account for the 
value of T c as well as the shape of the observed specific heat anomaly. By applying magnetic fields 
up to 6 T the hyperfine contributions Cm to the specific heat arising from the 55 Mn nuclei could be 
detected. From the time-dependence of the measured Cm the nuclear-spin lattice relaxation time 
Tin could be determined for the same field range in the temperature region 0.2 < T < 0.6 K. The 
nuclear magnetic relaxation was further studied by high field 55 Mn pulse NMR measurements of 
both the nuclear Ti n and T2 n at T — 0.9 K (up to 7 T). The data are in good mutual agreement 
and can be well described by the theory for magnetic relaxation in highly polarized paramagnetic 
crystals and for dynamic nuclear polarization, which we extensively review. The experiments provide 
an interesting comparison with the recently investigated nuclear spin dynamics in the anisotropic 
single molecule magnet Mni2-ac. 

PACS numbers: 75.10.Jm, 75.30.Kz, 76.60.-k 



I. INTRODUCTION 



Experimental and theoretical interest in magnetic 
molecular clusters carrying a net high spin has rapidly 
e volv ed in recent years (for reviews see, e.g., Refs. 
I ll l'2l.{l l lllfil) . Since the cores of these molecules can be 
viewed as nano sized pieces of magnetic insulators, they 
offer attractive possibilities to study magnetic objects of 
dimensions in between atom and bulk. Of great impor- 
tance is the fact that these (macro)molecules form stoi- 
chiometric chemical compounds, which may crystallize as 
molecular crystals, implying that for a given compound 
identical magnetic molecules arc arranged on the sites 
of a regular three-dimensional lattice. More often than 
not, there is only a single molecular site per unit cell, so 
that the symmetry axes of all molecules in the crystal 
are perfectly aligned. Provided that the intermolecular 
magnetic interactions are sufficiently weak, macroscopic 



solid state techniques can then be exploited to study the 
properties of individual cluster spins, taking into account 
their couplings to the "environment" (phonons, nuclear 
moments) as perturbations. Accordingly, such experi- 
ments have already provided highly interesting informa- 
tion about the quantum tunneling properties of the clus- 
ter spins 
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On the other hand, the long-range magnetic ordering 
(LRMO) phenomena expected to occur at sufficiently low 
temperatures as a consequence of the intercluster mag- 
netic dipolar coupling present an interesting object of 
study in itself^ii^ In many of these compounds, the clus- 
ters are only or mainly coupled by van der Waals forces in 
the molecular crystal. Short-ranged superexchange inter- 
actions may then be neglected, leaving only the dipolar 
coupling between cluster spins as a source for producing 
LRMO. Since in the literature of magnetic phase transi- 
tions few examples of LRMO produced purely by dipo- 



lar forces are yet available ji2i22i2i the chance to exploit 
molecular magnets to this end is quite attractive and 
could represent an important contribution to this field. 

However, for most high-spin molecules studied so 
far, such as Mn 12 ^ Feg^ 3 ., and the Mn 4 family of 
compoundst 2 ^ 2 ^^ the cluster spins have strong Ising- 
type anisotropy, associated with the zero-field splitting 
(ZFS) of the magnetic energy levels by the action of the 
crystal field. As a consequence, the cluster spins become 
frozen below a blocking temperature, Tb, of typically a 
few K, with the spin direction randomly distributed be- 
tween the two possible orientations along the easy axis. 
Evidently, this superparamagnetic blocking process is in 
competition with intercluster magnetic interactions that 
tend to establish LRMO at (usually much) lower temper- 
atures. Although it has been shown theoretically 2 ^ that 
the occurrence of quantum tunneling between opposite 
spin directions at temperatures below Tb may in prin- 
ciple produce sufficient fluctuations to overcome block- 
ing, in most of the investigated anisotropic molecules the 
times involved for the actual observation of the ensuing 
LRMO are still much too long. The recently discovered 
example of LRMO found at T c = 0.21 K in the molecular 
magnet Mn 4 Me is an exception rather than the rulei 2 ^ 

The obvious route to find dipolar-induced LRMO in 
molecular magnetic cluster compounds is, therefore, to 
search for high-spin molecules with as low anisotropy 
as possible and with negligible superexchange interac- 
tions. In a preliminary report 2 ^ on the compound 
Mne04Br4(Et2dbm)6 (hereafter called Mn§) we could 
show that it provides an excellent example. The Mn6 
molecule has a highly symmetric cluster core, comprising 
an octahedron of Mn 3+ ions the faces of which are capped 
by O 2- or Br~ ions. The structure of the molecule^ and 
a sketch of its octahedral core are shown in Fig. ^ From 
previous magnetic studies^ above 2 K it was found that 
the superexchange paths formed between the Mn 3+ ions 
(each having atomic spin s = 2) through the intervening 
2 ~ and Br~ ligands result in a relatively strong ferro- 
magnetic interaction, of value Jf/k-Q ~ +13 K on basis of 
the pair Hamiltonian TL = — 2JSi ■ Sj. As a consequence, 
the ground state is a S = 12 multiplet and the energy of 
the nearest excited state is approximately 150 K higher. 
The unit cell is monoclinic, with space group Pc, and 
contains four molecules that have such a high (nearly Td) 
symmetry that the net anisotropy for the cluster spin is 
quite small. No superexchange paths connecting neigh- 
boring clusters can be discerned indeed in the crystal 
structure, so that we can safely assume that the crys- 
tal binding arises solely from Van der Waals bonds. We 
note that, although intercluster magnetic ordering has 
also been reported for the molecular magnets Feig^i and 
Mn 4 Br, 32 in those cases superexchange between clusters 
apparently plays an important role, as evidenced, e.g., by 
the much higher T c values found (1.2 — 1.3 K). However, 
for a Cr 4 S cluster^ for which the intracluster exchange 
between the Cr 3+ ions happens to be likewise ferromag- 
netic (net spin 5 = 6), the low value of T c = 0.17 K that 
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FIG. 1: (a) Structure of Mn 6 04Br 4 (Et2dbm)6; (b) sketch 
of the symmetric octahedral core, containing six ferromag- 
netically coupled Mn 3+ ions, yielding a total spin S = 12 
for this molecular superparamagnetic particle. The gray and 
black spheres indicate the Br - and 2 ~ ions respectively. 
The white arrows illustrate the ferromagnetic alignment of 
the Mn 3+ s — 2 spins inside the cluster. 

was observed could be compatible with dipolar-induced 
magnetic ordering, in this case of antiferromagnetic type. 
More data would be needed, however, to substantiate 
this. 

In a preliminary report^ we could already show that 
the magnetic anisotropy in Mn@ is sufficiently low to en- 
able measurements of its magnetic susceptibility and spe- 
cific heat under thermal equilibrium conditions down to 
the lowest temperatures reached (T ~ 15 mK) by our 
experimental setups. The data did evidence a transition 
to LRMO, as hoped for, at a temperature of T c = 0.16 
K, corresponding apparently to a ferromagnetic arrange- 
ment of the cluster moments. Comparison with Monte 
Carlo simulations strongly supported the expected dipo- 
lar origin of the intercluster coupling. When applying 
magnetic fields of up to 6 T, the study of the time- 
dependent magnetic specific heat revealed a transition 
from fast relaxation to non equilibrium conditions within 
the experimental time window of 1 — 100 s. 

In the present work we extend these experimental and 
theoretical studies and discuss in much more detail the 
results and conclusions. In addition, we have performed 
55 Mn NMR studies in varying field, enabling to draw 
more definite conclusions about the magnetic relaxation 
of both electronic and nuclear spin systems in this ma- 
terial. The nuclear spin-lattice relaxation time is gov- 
erned by fluctuations of the cluster electronic spins, and 
is indeed quite fast in zero field. By applying a mag- 
netic field these fluctuations become progressively sup- 
pressed as a consequence of the Zeeman splitting of the 
electronic energy levels, thereby bringing the nuclear 
spin system out of thermal equilibrium. This provides 
an interesting comparison with recent zero-field Mn 
NMR studies 3 ^* 3 ^ 3 ^ 3 ^ of the anisotropic molecular mag- 
net Mni2-ac, for which the suppression of the magnetic 
relaxation in the thermally activated regime can be fully 
ascribed to the strong splitting of the cluster spin levels 
by the crystal field. For the present compound, crystal- 
field splittings play a very minor role in the relaxation 
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process, except for providing the necessary channel for 
energy transfer between spins and lattice phonons. We 
present a full analysis of the longitudinal and transverse 
nuclear relaxation in terms of previously developed the- 
ories for relaxation by paramagnetic impurities and for 
dynamic polarization, taking into account electron spin 
fluctuations by both spin-lattice relaxation and spin-spin 
interactions. Our NMR data in high fields prove to be 
in excellent qualitative as well as quantitative agreement 
with such theoretical predictions. The values for the lon- 
gitudinal nuclear spin-lattice relaxation rate and for the 
effective hyperfine interaction constant deduced from our 
high-field time-dependent specific heat data are likewise 
in good accord with the NMR results. 

The outline of this paper is as follows. After giv- 
ing a few experimental details in the next section, the 
measured susceptibility and zero-field specific heat data 
are presented and discussed in section III A and B, fol- 
lowed by a Monte Carlo simulation study of the zero-field 
specific heat in section III C. In section III D the field- 
dependent specific heat measurements are discussed, fol- 
lowed by the nuclear resonance and nuclear relaxation 
data in section III E. Section IV contains an analysis of 
these data in terms of existing theoretical models for re- 
laxation in highly polarized magnetic systems. Conclud- 
ing remarks are given in section V, while in Appendix A 
we describe the calculation of the demagnetizing factor 
for powder samples. Systeme International, units will be 
used throughout the paper. 



sample inside one of the two secondary coils and feeding 
the primary with an ac current (typically ~ 50 /iA), the 
induced voltage across the secondary is proportional to 
the susceptibility of the sample. By phase-sensitive de- 
tection we can also discriminate the real and imaginary 
parts of the susceptibility. The excitation frequency v 
was varied between 230 and 7700 Hz. Additional mea- 
surements above 1.8 K were performed using the ac op- 
tion of a commercial Superconducting Quantum Interfer- 
ence Device (SQUID) magnetometer. 

As for the NMR experiments, we introduced the sam- 
ple and the four-turns NMR coil inside the plastic 3 He 
pot of a pumped 3 He cryostat, where variable magnetic 
fields up to 8 T could be applied by a superconducting 
magnet. The resonance of the 55 Mn nuclei of Mng was 
observed by means of the spin-echo technique, with a 
typical duration of the ir/2 pulse of ~ 2 /is. 

Estimates of the magnitude of the magnetic anisotropy 
have been obtained as follows. Magnetic data measured 
above 1.8 K (Ref. H| can be fitted by expressions valid 
for fully isotropic S = 12 spins, i.e. they do not evi- 
dence any detectable ZFS for the total spin. The data 
are excellently fitted by the Brillouin curve calculated for 
an isotropic paramagnet with S = 12 and g = 1.98 [as 
derived from electron paramagnetic resonance (EPR)]. If 
we write a single-spin Hamiltonian for the molecule as: 

H = -£S z 2 - 5M BB a -S, (1) 



II. EXPERIMENTAL DETAILS 

Polycrystalline samples of Mng were prepared as re- 
ported in Ref. I30l Low-temperature specific heat mea- 
surements were performed in a home made calorimeter 
that makes use of the thermal relaxation method. 5,15 - 38 
For the measurements, a few milligrams of sample were 
mixed with Apiezon N grease and placed on the sapphire 
plate of the calorimeter. Details of this measurement 
technique are given in Ref. 1381 An important advantage 
of this method is that the characteristic time r of the 
experiment (typically, t c ~ 0.1 — 1000 seconds at low T) 
can be varied (within limits) by changing the dimensions 
(and therefore the thermal resistance) of the Au wire that 
acts as a thermal link between the calorimeter and the 
mixing chamber of the dilution refrigerator. Magnetic 
fields up to 16 T can be applied with a superconducting 
magnet and the lowest temperatures reached are of the 
order of 50 mK. 

The ac-susceptibility measurements were performed 
between 15 mK and 4 K in a home made susceptometer, 
placed inside the plastic mixing chamber of a dilution 
refrigerator and thermalized by the 3 He flow,— The sus- 
ceptometer, based on the mutual inductance technique, 
consists of a primary coil with 250 turns of 0100 /um 
NbTi wire, and two oppositely wound secondary coils, 
each with 660 turns of 040 /im Cu wire. By placing the 



these experiments provide an upper limit of \D\/ks <■ 
0.01 KiSi In order to obtain an independent estimation of 
D, high-frequency EPR experiments were carried out at 
the NHMFL in Tallahassee by J. Krzystek, using several 
frequencies in the range 95 - 380 GHz. Simulations of 
the spectra performed using Eq. agree well with the 
experimental results for 0.01 K < \D\/k B < 0.05 K. A 
value of \D\/ks — 0.01 K seems therefore appropriate 
to describe the ZFS in Mng. When discussing the ac- 
susceptibility and specific heat data in sections IlII Al and 
IIII Bl we shall adopt the value D/k-Q = 0.013 K, which 
yields the best agreement between theory and experiment 
for both techniques. 

The isotropic character of the molecular spin might 
seem paradoxical at first, considering that the individ- 
ual Mn 3+ ions, being Jahn- Teller ions, experience strong 
anisotropy: typical values of \D\/kB for the ion are a few 
tenths of Kelvin. However, the net D value entering in 
the spin Hamiltonian for the cluster can be seen in first 
approximation to result from the vectorial addition of 
the local anisotropy tensors of the individual ions, which 
then can give rise to a low net anisotropy for highly sym- 
metric molecules such as Mn 6 (cf. Fig. Q), even for large 
ZFS of the constituting atoms^ In fact, the possibility 
to tune the net anisotropy of the cluster spin by means 
of molecular synthesis is one of the attractive properties 
of these nanosized molecular superparamagnets. 
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III. EXPERIMENTAL RESULTS AND 
ANALYSES 

A. Magnetic susceptibility 

Strong evidence for the long-range ordering of the mag- 
netic moments is provided by the magnetic susceptibility 
data, shown in Fig. The real part x' °f the complex 
ac-susceptibility is plotted in Fig. |2Ja), and is seen to 
show a sharp maximum at T c = 0.161(2) K. We first 
demonstrate that the value of x' at T c is of the order of 
the estimated limit 1/N e g for a ferromagnetic powdered 
sample, where N e g is an effective demagnetizing factor 
appropriate for the (cylindrically shaped) container filled 
with the grains. In the Appendix we argue that N e g can 
be approximately given in terms of the demagnetizing 
factors iVg ra i n and iV con t of, respectively, the individual 
grains and the container as [Eq. l|A3(l ]: 

iVeff = iVgrain + /(JVoont ~ 1/3), (2) 

where / denotes the volume filling fraction of the con- 
tainer. 

Assuming the shape of the grains to be approxi- 
mately spherical, we put iV gra i n = 1/3, while from the 
shape of the container, we estimate N con t — 0.2. The 
density of the material is estimated to be p gra in — 
1.45 g/cm 3 from the value for the similar compound 
Mng04Br4(Mc2dbm)6, and the filling fraction is esti- 
mated as / ~ 1/3. All this then leads to iV cff ~ 0.29(5) 
and therefore x'(T c ) ~ 1/N e g = 3.5(6), the large error 
arising obviously from all the uncertainties in the above 
line of argument and in the estimates of the parameters 
involved. Next we should realize that this value would 
be valid for the x' measured along the easy axis, whereas 
even a relatively small anisotropy will lower appreciably 
the x' along the other directions. 40 Therefore, the pow- 
der x' could easily be lower by a factor of 2 [see Eqs. 
l|A5)l and jX6)| of the Appendix]. 

Given all the uncertainties, the above derived value can 
obviously serve as an order of magnitude estimate only, 
but we note that the experimental value of x'(^c) — 3 is 
indeed rather close. As a second argument for the fer- 
romagnetic nature of the transition we include therefore 
in Fig. Ufa) the powder susceptibility expected for the 
paramagnetic (non interacting) calculated from 

the spin hamiltonian (Eq. ^) and applying in addition 
the corrections for demagnetizing effects as described in 
the Appendix. Besides the curve for D/k-Q = 0.013 K 
appropriate for Mng , we also include for comparison the 
fully isotropic (D — 0) and infinite anisotropy limits. 
From this plot it is evident that, when approaching T c , 
the susceptibility of Mng increases appreciably above the 
paramagnetic limit, confirming the ferromagnetic nature 
of the correlations. 

Below T c , the powder x' decreases rapidly, as expected 
for an anisotropic ferromagnet in which the domain-wall 
motions become progressively pinned. The associated 
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FIG. 2: (Color online) Real (a) and imaginary (b) compo- 
nents of the ac-susceptibility at the indicated frequencies. In- 
set: magnification of x'(T) to evidence the frequency depen- 
dence of the peak. The solid lines in panel (a) give the calcu- 
lated behavior of the paramagnetic susceptibility (no interac- 
tions) for three values of the crystal field anisotropy constant 
D. In these calculations, the effect of the demagnetizing fields 
of the grains and the sample holder has been introduced as 
described in the Appendix. 



domain-wall losses should then lead to a frequency de- 
pendent maximum around T c in the imaginary part, x", 
as the experimental data of Fig. EJb) indeed show. In 
fact, although the Mng spins can be considered as nearly 
isotropic at high temperatures, the anisotropy energy 
is still large compared with the dipolar interaction en- 
ergy /lofj, 2 /AirkBr 3 ~ 0.1 K between nearest neighbor 
molecules. Thus the ordering should be that of an Ising 
dipolar ferromagnet. 

As shown in detail in the inset of Fig. |2Ja), the temper- 
ature Tpoak at which the maximum value of x' is found 
depends only weakly on v, which we attribute to the 
anisotropy. The total activation energy for the reversal 
of each Mng molecular spin amounts to DS 2 ~ 1.9 K, 
i.e., about 35 times smaller than for Mni2- Although 
this is quite small, one could still expect the super- 
paramagnetic blocking of the Mng spins to occur when 
T ~ T B (Mni 2 )/35, that is below ~ 0.15 K. Since this 
value is very close to the actual T c , one may expect that 
for T — > T c the approach to equilibrium begins to be hin- 
dered by the anisotropy of the individual molecular spins. 
We stress, however, that the frequency dependence of x' 
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observed here is quite different from that of the well- 
known anisotropic superparamagnetic clusters. A way to 
quantify the frequency dependence of the peak in yf is by 
means of the parameter AT pca k/[T , pca ] < A(log 1 Q z/)], which 
gives the variation of T pea k per decade of frequency. We 
find here AT pcak /[T pcak A(log 10 v)] ~ 0.03 - 0.05, to be 
compared with the typical values of ~ 0.20 for super- 
paramagnetic blocking. In fact it is closer to the value 
~ 0.06 found for certain types of spin glasses^ but the 
peak observed here is much higher and sharper. Also, 
since the cluster spins are situated on a regular crys- 
tal lattice, a comparison with random magnetic systems 
would not be appropriate. The frequency dependence we 
observe is indeed very much weaker than is found in the 
LiHo 2 ,Yi_ :c F4 system with x = 0.045 (Ref. 42), for which 
x value that material is in the "antiglass" regime42*i4 In 
Mng we found that T vea ^(u) ~ T pca k(0) + Kv a , with 
a ~ 0.4 and a zero-frequency limit of the peak in x' 
Tpeak(O) — 158 niK. At essentially the same temperature, 
we find a fairly sharp peak in the zero-field magnetic spe- 
cific heat (see next section), instead of a broad anomaly 
as observed in LiHoo.o45Yo.95sF4. Also this finding ap- 
pears to exclude an interpretation in terms of a freezing 
transition in Mng. 

We finally turn to the susceptibility as measured above 
T c in the paramagnetic region, which was plotted as 1 /%' 
vs T in the inset of Fig. 2 in Ref. 2!). We first note that 
no evidence for relaxation effects were found in this range. 
Up to a frequency of 7700 Hz no appreciable x" was 
detected and the measured x' smoothly joins the data 
measured above 2 K with the SQUID susceptometer. We 
may therefore conclude that in the whole temperature 
range down to T c the spin-lattice relaxation time is quite 
short, of the order 10~ 4 s or less. 

The high-temperature susceptibility data x[ corrected 
for the demagnetizing field (x[ — x'/fi ~ A^ffX'D fol- 
low the Curie- Weiss law x[ — C/(T — 8) quite well 
down to approximately 0.3 K, with C = 0.62(2) K and 
9 = 0.14(3) K. The constant C equals, within the exper- 
imental errors, the theoretical value for isotropic spins 
N A g 2 ^%S(S + l)/3k B V m = 0.595 K, where S = 12, 
g = 2, and V m — 1647 cm 3 /mole is the molar volume. 
The positive 9 confirms the ferromagnetic nature of the 
ordered phase. The fact that mean field theory is so well 
obeyed down to very close to T c is as expected for a dipo- 
lar ferromagnetiaSiaSiaM^ We remark that the behavior of 
the powder susceptibility we observe for Mn 6 in Fig. |2 
closely resembles previous powder data for Cs2NaGdCl6 
(cf. Fig.3 of©, and LiHo^Yi-^ (x = 0.46) (cf. Fig.l 
in and Fig. 4 in 142^ . Both materials are considered 
to be examples of anisotropic dipolar ferromagnets. Al- 
though weak ferromagnetism (i.e., canted antiferromag- 
netism) would also lead to x(T = T c ) = 1/N, this would 
be accompanied by a negative Curie- Weiss 6, whereas we 
observe a positive value. Moreover, for canted antiferro- 
magnetism an antisymmetric interaction term of the form 
d • (Si x Sj) is needed (see e.g. Ref. l50|l and this is not 
expected for a system of equivalent magnetic moments 
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FIG. 3: (Color online) Zero-field specific heat of Mn6 plotted 
as a function of temperature. Solid line: phonon contribution 
{C/R = 0.010 T 3 ). Dashed-dotted line: Schottky contribu- 
tion due to crystal field splitting of the S — 12 multiplet as 
calculated with D/ks = 0.013 K. Dashed curve: expected 
nuclear contribution from the Mn nuclear spins. 



interacting by dipolar interactions, as is the case here. 
For metamagnetic systems, having ferromagnetic nearest 
neighbor interactions and weaker further neighbor cou- 
plings, also a positive value for 9 can be found but the 
magnetic ordering below T c due to the further neighbor 
interaction is basically antiferromagnetic. Consequently, 
although the value of %(T = T c ) can be much higher than 
for nearest neighbor antiferromagnets, it will fall still far 
below the ferromagnetic limit of 1/N. We therefore con- 
clude that the evidence for ferromagnetic dipolar order in 
Mn6 presented here is quite strong. Additional studies, 
for instance of the spontaneous magnetization below T c , 
would of course be welcome to provide further proof. 



B. Zero-field specific heat 

Specific heat data c taken in zero applied field are 
shown in Fig. [3] as a function of temperature on a dou- 
ble logarithmic scale. Above 2 K, the specific heat is 
dominated by the contribution from the lattice phonons 
that can be reasonably fitted by the well-known low-T 
Debye approximation: q a tt oc (T/8d) 3 , with a Debye 
temperature Qd — 29 K. Such low values are commonly 
observed in the molecular cluster compounds^4i^& re- 
flecting the weak bonding between the cluster molecules 
in such molecular solids. 

In the lowest temperature range below 0.1 K, the spe- 
cific heat is seen to remain rather high valued, which can 
be mainly ascribed to the expected contribution from the 
55 Mn nuclear spins, the energy levels of which are split 
by the hyperfine interaction with the Mn 3+ electronic 
spins s = 2 (see section III D below). This contribu- 
tion will become more clearly and directly visible in the 
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field-dependent studies, where it will be shown to be de- 
scribable by a term c nuc \T 2 /R ~ 5.2 x 10~ 3 K 2 , which 
is plotted as the dashed line in Fig. 03 Anticipating 
the discussion below, we subtract this hyperfine contri- 
bution, as well as the phonon oc T 3 term from the raw 
data in Fig. in order to obtain the magnetic specific 
heat c e i associated with the cluster electron spin S = 12 
only, as plotted in Fig. 01 The resulting curve shows two 
characteristic features. At a temperature T c ~ 0.15(2) 
K a peak is observed, that can be associated with the 
transition to long-range magnetic order, in good agree- 
ment with the value for T c deduced from the peak in the 
ac-susccptibility extrapolated to zero frequency. Above 
T c , one observes a widely extended "high-temperature 
tail" , that reflects the weak zero-field splittings (ZFS) of 
the S = 12 multiplets by the crystal field interactions. 
We recall that even for a D value as small as 0.013 K 
the total splitting of an S = 12 multiplet will still be 
an appreciable DS 2 — 1.9 K. In the absence of mag- 
netic intercluster interactions the ZFS of such a multiplet 
would lead to a multilevel Schottky curve, shown as the 
dashed-dotted curve in both Fig. [31 and 01 where the fit 
to the experiment leads to the estimate |D|//cb — 0.013 
K, in reasonably good agreement with the values quoted 
above. In contrast with the highly anisotropic molecular 
clusters Mni2, Mn4, and Fes, where the multilevel ZFS 
Schottky is found above 1 K and is the most pronounced 
feature of the experimentally observed magnetic specific 
heat; 5 i 15 i 28 its presence in Mng is masked at the low-T 
side by the ordering anomaly produced by the effects of 
the intercluster magnetic interactions. Numerical inte- 
gration of Cei/T between 0.08 K and 4 K gives a total 
entropy change As e i/R = 3.5(2) per mole, close to the 
expected total entropy for a fully split S — 12 multiplet, 
namely As cl /R = ]n(2S + 1) = 3.22. This confirms the 
consistency of the above subtraction procedure used to 
obtain c c \. The variation with temperature of the en- 
tropy, s e i/R, is also shown in Fig. 01 We note in partic- 
ular that at T c itself, the entropy only amounts to about 
1R per mole, indicating that only the lowest energy levels 
of the cluster spins are involved in the actual magnetic 
ordering process, the majority of the higher-lying lev- 
els being already depopulated (the lowest lying doublet 
ground state on its own would already give a contribu- 
tion to the entropy of i?ln2 = 0.69i? per mole). Indeed, 
the distance to the nearest lying excited state would be 
equal to (2S — 1)D ~ 0.3 K in the non interacting limit. 



C. Monte Carlo simulations 

To simulate the zero-field specific heat data, Monte 
Carlo (MC) calculations were performed for an S = 12 
Ising model of magnetic dipoles on a body-centered or- 
thorhombic lattice containing Z = 2 molecules, with 
axes a x = 15.7 A, a y = 23.33 A, and a z — 16.7 
A. This choice approximates the crystal structure of 
Mn604Cl4(Et2dbm)6, that is the most closely related 
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FIG. 4: (Color online) (a) T dependence of the electronic 
specific heat at zero applied field (circles) . The dashed-dotted 
line is the Schottky anomaly calculated with D/Icb = 0.013 
K. The solid line is the Monte Carlo (MC) calculation for an 
orthorhombic lattice of 1024 Ising spins with periodic bound- 
ary conditions. For each point, 2 x 10 4 MC steps per spin 
were performed, (b) T dependence of the entropy obtained 
by integration of the electronic specific heat curve. 



compound for which the structure could be resolved. To 
further simplify the Monte Carlo simulations, we approxi- 
mated the monoclinic structure by an orthorhombic one. 
The density, p = 1.45 g/cm 3 , was estimated from the 
Mn604Br4(Me2dbm)6 compound. The resulting molar 
volume is a few percent larger than what would be ob- 
tained from the lattice parameters used in the Monte 
Carlo simulation, but this discrepancy is due entirely 
to the approximation of orthorhombic structure. The 
Hamiltonian includes the dipolar interaction term as well 
as the anisotropy term —DS 2 given in Eq. (JIJ. For 
T < 0.5 K the intermolecular dipolar interactions be- 
come important and remove the degeneracy of the | ± m) 
spin doublets. The MC simulations show that the ground 
state is indeed ferromagnetically ordered, as observed, 
and predict a shape for c e i that is in good agreement with 
the experiment. In the upper panel of Fig. 01 we show 
c i as calculated assuming all molecular easy (z) axes to 
point along a z , i.e., one of the two nearly equivalent short 
axes of the actual lattice. Similar results were obtained 
for other orientations chosen for the anisotropy (z) axis. 
We note that our simulations give T c = 0.22 K, which is 
slightly higher than the experimental T c = 0.161(2) K. 
This difference may be due to the Ising approximation 
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FIG. 5: Calculated critical temperature T c as a function of 
the anisotropy parameter D, for classical Heisenberg spins. 



taken for the intercluster dipolar interaction and to un- 
certainties in the values of the lattice parameters. In fact 
we note in passing that an almost perfect coincidence of 
our calculated curve with the experimental data may be 
obtained by assuming a smaller value of the magnetic 
moment, i.e., taking g = 1.75 instead of g = 2.00. 

To pursue this point further, additional MC simula- 
tions were performed for the same crystal lattice, but 
now with classical Heisenberg spins replacing the S = 12 
Ising spins. To investigate the sensitivity of our results 
to the type of anisotropy, the sign and magnitude of D 
were varied. These calculations resulted in the phase di- 
agram shown in Fig. [3 which we include here since it 
illustrates the complicated way in which the nature of the 
actual ground state and the value of T c may depend on 
the combination of the long-range dipolar interaction and 
the anisotropy parameter. Although the ground state for 
this lattice is always found ferromagnetic, it can be ei- 
ther "uniaxial" , with strong preference for the spins to lie 
along the a y axis [chosen for this particular example to 
be the anisotropy axis z of Eq. Q], or "planar", in the 
sense that the a x — a z plane becomes an easy plane, with 
a weak preference for a given direction in the plane, as 
sketched in the figure. Interestingly, the switching point 
between these "uniaxial" and "planar" orientations is not 
at D = 0. The reason for this is as follows. Because the 
crystal lattice under consideration is far from being cu- 
bic, the dipolar interaction energy is rather anisotropic. 
The dipolar energy is minimized when the magnetization 
M points in the direction (in the a x — a z plane) shown in 
the inset on the left-hand side of Fig. [3] Therefore, M 
points along this direction for D < 0. On the other hand, 
the energy minimization for D > is a competing pro- 
cess. Clearly, dipolar interaction must become dominant 
for sufficiently small values of D. The numerical results 
show that this occurs if < D/k-Q < 3 mK. The numeri- 
cal datapoints in Fig. [SJalso show that T c varies sharply 



within the —0.01 < D/k-Q < 0.01 K range. System size 
effects and computer time restrictions do not allow us to 
determine whether T c vanishes completely. The lowest 
numerical value obtained is as small as T c ~ 0.03 K at 
D/fcfl — 3 mK. Outside this narrow range of D, T c is 
already almost equal to the limiting values of ~ 0.7 K 
and ~ 0.3 K, reached for infinite negative and positive 
D, respectively. Such a variation of T c with anisotropy, 
as well as the form of the calculated and observed specific 
heat ordering anomaly, appear to be specific for dipolar 
interactions. They differ widely from the correspond- 
ing behavior known for three-dimensional (Heisenberg, 
Ising, XY) ferromagnetic lattices with nearest-neighbor 
interaction only^i For instance, in those models the vari- 
ation of T c with anisotropy is restricted to about 20% 
of the T c value, which moreover is highest for the Ising 
("uniaxial") case. It is also interesting to compare these 
Monte Carlo calculations with the predictions of simple 
mean field theory. The latter gives T c = 2S 2 J c g/kB and 
T c = 25(5+ l)J c ff/3fcB, respectively for D — ► +oo and 
D = 0, where J c g is an effective interaction constant. 
The ratio between these two limits, which is 35/(5 + 1), 
is about 2.5 times larger than what is obtained by Monte 
Carlo calculations. 

In concluding this section we may stress that our de- 
tailed calculations have evidenced that for a dipolar mag- 
net the value of T c may vary strongly with anisotropy and 
lattice symmetry. This illustrates the danger in drawing 
conclusions about the nature of the magnetic interac- 
tions, i.e., whether they are of dipolar origin or not, just 
by comparing the value of koT c with the dipolar interac- 
tion energy of a pair of nearest neighboring spins. Our 
detailed calculations specific for this compound show a 
good quantitative agrement assuming just dipolar cou- 
pling, both as regards the value of T c and the shape of 
the specific heat anomaly. Furthermore, the prediction 
that the ordering is ferromagnetic, as observed, is quite 
robust since it is independent of the details of the simu- 
lations. All this confirms that LRMO in Mn@ is mainly 
driven by dipolar interactions. 



D. Field-dependent specific heat: Nuclear spin 
contribution 



We next discuss the time-dependent magnetic specific 
heat c m = c — ci a tt measured under varying applied mag- 
netic fields B a , as plotted in Fig. [5] Even for the low- 
est _B a value, the ordering anomaly is already fully sup- 
pressed, as expected for a ferromagnet. Accordingly, we 
may account for these data with the Hamiltonian Q ne- 
glecting dipolar interactions. The Zeeman term splits 
the otherwise degenerate | ± m) doublets, and already 
for B a ~ 0.5 T the level splittings become predominantly 
determined by B ai so that the anisotropy term can then 
also be neglected. As seen in Fig. [5] the calculations 
performed with D = (dotted curves) reproduce quite 
satisfactorily the data at higher temperatures. 
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FIG. 6: Temperature dependence of the magnetic specific 
heat c m at B a = 0.5 T (•), 1 T (□), 2 T (a), 2.5 T (o), 
3 T (*), and 6 T (o). The dotted lines represent the calcu- 
lated electronic contribution, whereas the dashed line is the 
expected nuclear specific heat at zero field calculated with 
Eq. The solid lines represent the total time-dependent 

c m (electronic + nuclear) calculated accounting for the nu- 
clear Ti n (see text). Inset: detail of c m (T) at B a = 2 T, for 
long (~ 100 s) and short (~ 1 s) experimental times. The 
dashed lines are the calculated contributions arising from the 
55 Mn and 1 H nuclear spins. 



(see section III E), we obtain the dashed line shown in 
Fig. at zero field. This is the same contribution that 
was subtracted from the zero-field data shown in Fig. 01 
The difference between the calculated and experimental 
c nuci> wmc h becomes especially evident for the B a = 2T 
curve, can be due to the shift of the 55 Mn nuclear energy 
splittings by the applied field, which is neglected in Eq. 

A remarkable feature of the experimental data that is 
not reproduced by these equilibrium calculations is that, 
at the lowest T, the nuclear specific heat drops abruptly 
to a baseline of about 3 x 10 _3 i?. This remaining specific 
heat is probably a background feature arising from in- 
complete correction for the field-dependent addenda con- 
tributions. The crossover temperature T* where the drop 
of c m occurs depends on B a but also on the characteristic 
time constant r e of our (time-dependent) specific heat ex- 
periment: as is shown in the inset of Fig. [f)|the deviation 
from the (calculated) equilibrium specific heat is found 
at a lower T when the system is given more time to re- 
lax. Interestingly, the specific heat becomes even smaller 
than the expected contribution of the protons. The drop 
therefore shows that, below T*, the nuclear spins of both 
the 1 H and 55 Mn atoms cannot attain thermal equilib- 
rium with the lattice phonons within the experimental 
time r e because the longitudinal nuclear spin-lattice re- 
laxation time becomes too short. This relaxation effect 
can be described as follows: 



However, when the maxima of the multilevel Schot- 
tky anomalies are shifted to higher T by increasing B a , 
an additional contribution is revealed at low T. It is 
most clearly visible in the curves for 1 T < B < 2.5 
T, and varies with temperature as c m T 2 /R ~ 4 x 10~ 3 . 
We attribute this contribution to the already-mentioned 
high-temperature tail of the equilibrium nuclear specific 
heat c^~]. As discussed in section IV, this specific heat 
should be dominated by the contributions of the six 55 Mn 
nuclear spins (/ = 5/2) and of the 114 protons present 
in each molecule. The energy levels of the former are 
split, even at zero field, by the strong on-site hyperfine 
interactions with the Mn 3+ electronic spins s = 2. For 
the specific heat analysis this interaction can be approx- 
imated by Hm — AI ■ s, where A is an effective isotropic 
hyperfine coupling constant. By contrast, we may neglect 
the hyperfine splitting of the protons because it can be 
expected to be small compared to their nuclear Zeeman 
splitting for B a > 1 T. This is indeed confirmed by the 
NMR experiments shown below. The high-temperature 
limit of the nuclear specific heat can therefore be approx- 
imated by the sum of two contributions: 52 



„(eq)r 



'nucl^ 



R 



6 x 



A 2 s 2 I(I+ 1) + 114 x 



2fc B 



(3) 



where 7h = 2.675 x 10 s rad T _1 s _1 is the protons gy- 
romagnetic ratio. Taking A/k^ = 8.6 mK as estimated 
from the NMR spectra measured for the same sample 



^) = Cl ( 1 C u q c ) Jl-cxp(-r c /T ln )], 



(4) 



where, to simplify the discussion, we have used the same 
nuclear spin-lattice relaxation (NSLR) time Ti n for both 
protons and 55 Mn although they can obviously differ 
from each other. According to Eq. J3J the specific heat 
decreases fast when Ti n becomes of the order of r e . This 
crossover to a non equilibrium regime (as measured by 
time-dependent specific heat) provides therefore direct 
information on the temperature and field dependence of 
the nuclear T ln . 

As is well established^ Ti n (T,B a ) can be related to 
the time-dependent part of the transverse hyperfine field 
as produced by the fluctuations of the electron spin. For 
the case of Mn 6 both the electron spin fluctuations due to 
spin- lattice coupling and to dipolar spin-spin interactions 
will have to be considered. A more extensive theoretical 
treatment is given below in section IV, in terms of ex- 
isting models for nuclear relaxation in magnetic crystals. 
As will be seen, this treatment predicts the behavior of 
the longitudinal nuclear spin-lattice relaxation time at 
low temperature and high fields (> 3 T) to be given by 



1 d3 ( 9mB z 

— « K0 B a exp f-_ 



(5) 



where kq is a constant that depends on the electronic 
spin-lattice relaxation rate and on the details of the re- 
laxation mechanism. This shows that an exponential 
temperature dependence of the NSLR rate is expected 
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FIG. 7: (Color online) Field dependence of the magnetic 
specific heat at T — 0.3 K (solid dots) and T = 0.6 K (open 
dots). The dashed and dotted lines represent respectively the 
calculated nuclear (see Eq. £J) and electronic contributions 
to the equilibrium specific heat. The thick and thin solid lines 
give respectively the equilibrium c m and the time-dependent 
c m calculated accounting for the field-dependent nuclear Ti n . 
Inset: B*(T*) obtained from T sweeps, as in Fig. |S| (open 
dots), or from B sweeps, as in the present figure (solid dots). 
From the fit (solid line) to Eq. we extract reo — 40 s _1 T~ 3 . 

at high fields. The effect of the field is to polarize the 
electron spins, which reduces the fluctuations of the hy- 
perfine field, thus effectively disconnecting the nuclear 
spins from the lattice. 

As seen from Eq. (JSJ the nuclear spins can be taken 
out of equilibrium either by decreasing T down to T* 
at constant field (as in Fig. |SJ|, or by increasing B a up 
to a given value B* at constant T. The latter effect is 
indeed also observed experimentally, as shown in Fig. [7| 
for T = 0.3 K and 0.6 K. In this figure the transition to 
non equilibrium is obvious from the fact that the data 
measured in high field fall far below the expected nuclear 
contributions (dashed lines). The inset of Fig. [7| shows 
B*(T*) obtained either from T sweeps at constant B (as 
in Fig. |BJ or from B sweeps at constant T (as in Fig. 
[7J|. The two methods prove to be fully consistent with 
each other. The fit of B*(T*) using Eq. (J5J gives an 



M(t) 
M(oo) 



where 1/Ti„ is the nuclear spin-lattice relaxation rate, 
and a is a stretching exponent, which needs to be in- 
troduced to account for the large inhomogeneity of the 
NMR line, which causes the inversion recovery to consist 



average value of kq ~ 40 s -1 T~ 3 . Using this value of 
kq we have calculated the time-dependent c nuc i from Eq. 
Q . Adding this to the calculated electronic specific heat 
yields the solid lines in Figs. [S]and[7| which can be seen to 
be in reasonably good agreement with the experimental 
data over the whole range of field and temperature. 



E. 55 Mn - NMR and nuclear relaxation 

The nuclear spin-lattice relaxation in Mng has been 
further investigated by 55 Mn nuclear magnetic resonance. 
Apart from the fundamental question as to by which 
mechanisms the magnetic relaxation proceeds in this 
isotropic molecular magnetic crystal, it is of interest to 
compare with the typical behavior observed recently for 
the highly anisotropic single-molecule magnets like Mni2- 

„ „ 34.35.36 

As is well known, it is difficult to observe NMR for 
nuclei of paramagnetic ions due to the very large and 
strongly fluctuating magnetic fields produced at the nu- 
clei by the electron spin through the (on-site) hyperfine 
interactions. As a consequence, nuclear resonance lines 
become very broad and spin-lattice relaxation rates too 
fast to be measured. To enable the observation of the 
NMR signals, one should therefore take recourse to the 
low-temperature and high-field regime, in which electron 
spin fluctuations can be expected to be sufficiently sup- 
pressed. Accordingly, we performed our experiments at 
T = 0.9 K, using a 3 He cryostat, and fields in the range 
3 to 7 T. 

The 55 Mn NSLR was studied by measuring the re- 
covery of the nuclear magnetization after an inversion 
pulse. By integrating the echo intensity we obtained 
recovery curves as those shown in Fig. OJa). For the 
ease of comparison between different curves, we renor- 
malize the vertical scale such that M(Q)/M(oo) = — 1 
and M(t > T ln )/M(oo) = 1, even though usually 
|M(0)| < |M(oo)|: this is just an artifact that occurs 
when the NMR line is much broader than the spectrum 
of the inversion pulse, and does not mean that the length 
of the 7T pulse is incorrect. Since the 55 Mn nuclei have 
spin / = 5/2, the recovery of the nuclear relaxation for 
the central line in the quadrupolar split manifold is de- 
scribed by^ 



(0) 

I 

of a combination of recoveries with different rates. We 
typically found an optimal value of a ~ 0.5, although 
the choice of the stretching exponent does not strongly 
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FIG. 8: (Color online) Inversion recovery (a) and decay of 
transverse magnetization (b) for the Mn nuclei at T = 0.9 
K, B a = 5 T and v = 251.5 MHz. The lines in (a) are fits to 
Eq. © with a ~ 0.5 (dashed) and a — 1 (solid). The solid 
line in (b) is a fit to Eq. 




FIG. 9: (Color online) 55 Mn NMR spectra at T = 0.9 K and 
measuring frequencies u = 246.5 MHz (down triangles) and 
v — 251.5 MHz (up triangles). The gap in the data around 
B a — 5.8 T is due to the crossing with the H line arising 
from the Zeeman-split proton levels. The lines are Gaussian 
fits with total width 2a B ^ 2.2 T. 



influence the value of 1/T ln extracted from the fit. 

The transverse spin-spin relaxation (TSSR) rate 1/T2 n 
was studied by measuring the decay of echo intensity 
upon increasing the waiting time r between the ir/2- 
and the 7r-pulses. The decay of transverse magnetization 
M±(t) can be fitted to a single exponential: 



M±(2t) 
M ± (0) 



,-2r/T 2n 



(7) 



as shown in Fig. [S^b). 

The field-sweep NMR spectra in Fig. clearly show 
that it is impossible to determine whether there are in- 
equivalent sites in the molecule, as regards the hyper- 
fine coupling [compare with the case of Mni2-ac (Ref. 
|{55|) ]. This may be due to the large quadrupolar split- 
ting expected in Mn 3+ sites, plus the fact that our sam- 
ple is an unoriented powder. As expected from the 
internal ferromagnetic structure of the cluster electron 



spins, the 55 Mn spectrum shifts to higher fields when 
lowering the frequency. The spectra can be fitted by a 
Gaussian shape with total width 2<tb — 2.2 T. If this 
width were due to quadrupolar splitting Aisq only, one 
would deduce Avq ~ 7 MHz: this estimate can be ob- 
tained from a comparison to the Mn' 1 ' line in Mni2-ac, 
where Av Q = 0.72 MHz yields 2a v = 2.4 MHz,S i.e. 
2<j b = 2ov/(7m„/27t) ~ 0.23 T, where 7 M „ = 6.64 x 10 7 
rad T _1 s _1 is the 55 Mn gyromagnetic ratio. Such an 
estimate is thus even larger than the highest Avq ~ 4.3 
MHz found in the Mn' 2 ' sites of the less symmetric Mni2- 
ac cluster S We expect therefore that the random orien- 
tation of the crystallites and, eventually, the presence 
of inequivalent Mn sites as regards the hyperfine cou- 
pling, are also contributing to the observed broadening. 
Indeed, when decreasing the frequency by 5 MHz the 
maximum of the spectrum shifts only by 0.24 T, instead 
of the 0.47 T that would be expected when all the local 
hyperfine fields are antiparallel to B a . We conclude that 
the observed spectrum, as well as the NSLR and TSSR 
data, should be considered as obtained from a mixture 
of nuclear signals arising from randomly oriented crystal- 
lites with largely overlapping and quadrupolar-split NMR 
lines from all the Mn sites in the cluster. Extrapolating 
to B a = the field dependence of the peak of the spec- 
trum, one obtains v(0) ~ 360 MHz =>■ i?h yP — 34 T, very 
similar to the highest value found in the Mn' 3 ' site of 
Mn 12 -aci££ Finally, in connection with the discussion of 
the nuclear specific heat of the previous section, it is in- 
teresting to mention that the *H resonance is found at a 
value of field given simply by £?h — w /7h (the excluded 
region in the spectra shown in Fig. EJ. This confirms that 
the local hyperfine fields do not appreciably shift the : H 
resonance frequency and therefore, for sufficiently high 
fields (-B a > 1 T), the hyperfine interaction of protons 
can be neglected, as we did. 

Figure ^| shows the field dependencies of the NSLR 
rate 1/T" ln and the TSSR rate 1/T 2n , measured at con- 
stant frequency v — 251.5 MHz and temperature T = 0.9 
K. From the discussion above it is clear that these data 
must be interpreted with a certain caution, since shifting 
_B a at constant v means that we are sampling each time 
a different portion of the NMR signal, which means dif- 
ferent quadrupolar satellites, different orientation of the 
crystallites, etc. Nevertheless, the agreement with the 
estimate of Ti„ obtained by specific heat data (inset of 
Fig. [fjj) turns out to be satisfactory. We can directly com- 
pare the NSLR rates 1/T in (i3 a ) obtained from respec- 
tively NMR and specific heat data by plotting (dashed 
line in Fig. 1/Ti n (£? a ) as calculated from Eq. ©, 
with T = 0.9 K and fixing k q = 40 s _1 T~ 3 as obtained 
from the fit of B*(T*) in the inset of Fig. The agree- 
ment is seen to be reasonable. The solid line in the same 
figure represents a fit to the l/Ti„(B a ) based on a model 
described by Eq. I|34l) of the following section, with pa- 
rameters given in the discussion there. 

In Fig. Uniwe also show the TSSR rate 1/T 2n (-B a ), with 
ordinate axis shifted in order to compare its field depen- 
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FIG. 10: (Color online) Field dependencies of the 55 Mn 
NSLR (full squares, left scale) and TSSR (open circles, right 
scale) at T = 0.9 K and v = 251.5 MHz. Notice the factor 
10 4 between left and right scales. The hatched area around 
5.8 T indicates the region where the 1 H line overlaps with the 
55 Mn resonance. Dashed line: calculated 1/Ti n (-B a ) accord- 
ing to Eq. Jjj} for the value ko = 40 s~ lr T~ 3 as extracted from 
specific heat experiments. This represents exactly the same 
line as in the inset of Fig. |7| Solid line: fit to Eq. 1M It below 
with parameters discussed in the text. Dotted line: 1/T2 n 
calculated according to Eq. 1251 . 



dence to that found for the NSLR rate. As explained 
below, we expect l/T2 n (S a ) cx exp^/iB^a/feeT) at high 
fields. This dependence, and the quantitative result cal- 
culated according to Eq. I|25[) , are shown as a dotted line 
in Fig. 1101 and yield the right order of magnitude of the 
observed T2„ ~ 0.1 — 1 ms. 



IV. DISCUSSION: MAGNETIC HYPERFINE 
INTERACTIONS AND NUCLEAR MAGNETIC 
RELAXATION 

In this section we apply theoretical results for magnetic 
hyperfine interactions and nuclear magnetic relaxation in 
magnetic insulating solids to interpret the data of the pre- 
vious sections. We recall that for such materials direct 
relaxation channels such as quadrupolar interactions con- 
necting nuclear spins to the lattice become ineffective at 
low temperatures. Thus the electron spins present have 
to serve as an intermediary between nuclei and phonons 
in some way or the other and we have to consider the in- 
terconnected spin dynamics of both nuclear and electron 
spin systems and their coupling to the lattice. As will be 
shown below, for highly polarized electron spin systems, 
such as Mng in high fields at low temperature, the relax- 
ation behavior can be adequately described in terms of 
theoretical models previously developed in the field of dy- 
namic nuclear polarization. These same models should 
also provide a good basis to describe the relaxation in 
highly anisotropic single-molecule magnets like Mni2-ac 
and Fes below their blocking temperatures, where an ex- 
treme polarization of the electron spins is induced by the 



crystal field. Although several groups have described the 
application to molecular magnets of some general the- 
ories of nuclear relaxation)2£i££i£LS it appears that the 
role of the electronic dipolar coupling has so far been 
overlooked. We therefore consider worthwhile to give a 
detailed overview of the different ingredients needed to 
arrive at a consistent picture describing the behavior we 
observe in both the NMR and the specific heat of Mng, 
being confident it will be quite useful for the other ma- 
terials mentioned as well. We will briefly come back to 
this point in the conclusion section, drawing a compari- 
son with our observations in Mni2-ac. We start with an 
evaluation of the relevant hyperfine interactions between 
the cluster spins and the various nuclear spins in the Mng 
molecular cluster. Quite generally^SS the hyperfine in- 
teraction hamiltonian of a nuclear moment I with the 
surrounding electron spins Sj can be written in the form 
of the bilinear coupling 

W hf = £>i- Sj (8) 

j 

Here the hyperfine interaction A is a second rank tensor 
and summation is over electron spins on both the same 
atom and surrounding atoms. It is often convenient to 
interpret the hyperfine coupling in terms of a magnetic 
hyperfine field: -Bhf = — (7n^) s • A, acting on the nu- 
clear spin in addition to the applied field B a . 7„ is the 
nuclear gyromagnetic ratio. The total magnetic Hamil- 
tonian for the nuclear spins then becomes 

Whf = -Inhl • (B a + B hf ) = - 7n « • B tot (9) 

In principle Hhf will be time-dependent since both A and 
sj can depend on time due to, respectively, atomic mo- 
tions (not considered here) and fluctuations of the elec- 
tron spins. Provided that the frequency of these fluctu- 
ations is fast compared to the nuclear Larmor frequen- 
cies produced in the static case, the nuclear resonance 
is still well defined be it at a frequency that is shifted 
with respect to that for Bhf = 0, the shift being pro- 
portional to the time average of of A{t) ■ s. In the ef- 
fective field picture the hyperfine field can be split up 
into a static part, (Bhf), and a time-dependent part 
Bhf(i) = (Bhf) + b(i). The time-dependent fraction, 
b(i), is usually much smaller and can then be treated as 
a perturbation that may produce relaxation of the nu- 
clear polarization. Neglecting the quadrupolar interac- 
tions, the remaining magnetic hyperfine interaction may 
be decomposed into contributions coming from the cou- 
pling of the nuclear moment with the orbital motions of 
the electrons, the dipolar interactions with the electron 
spins and, in case a finite density of electrons is present 
at the nuclear site, the part due to the Fermi-contact in- 
teraction. Thus we may write (Bhf) = Bdi P + Bl + Bp, 
in an obvious notation. As for the relative strengths of 
these contributions, for nuclei (such as the present 55 Mn) 
residing on magnetic atoms the Fermi-contact term is 
strongest by far, with B-p ~ 10 — 10 2 T, followed by 
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the orbital and dipolar interactions with electrons on the 
same atom (on site), typically of order 10 T. For ions 
with closed or half-filled shells both £?l and £?dip vanish, 
whereas also for non S-state ions, such as Mn 3+ , the or- 
bital moment can be quenched by crystal-field splittings, 
so that -Bl becomes negligible. Furthermore, nuclei on 
nonmagnetic ligand atoms directly coordinating the mag- 
netic ions via covalent bonds may also experience sub- 
stantial Fermi-contact interaction, of order 1 — 10 T. Nu- 
clei on more distant nonmagnetic atoms will mainly ex- 
perience the long-range dipolar interactions, giving typi- 
cally B dip - 0.1 - 1 T. 

For the 55 Mn nuclei on the Mn 3+ ions (with nuclear 
spin 7 = 5/2 and electron spin s = 2) of the present 
compound Mn6 we deduced in section III-E a resonance 
frequency (extrapolated to B a = 0) centered around 360 
MHz (corresponding to (-Bhf) = 34 T), which is almost 
the same as the value 364 MHz Kubo et alm^ recently 
found for one of the two Mn 3+ sites in Mn 12 -ac. Accord- 
ing to their analysis, the corresponding Bhf = 34.5 T for 
this site results from the combination of Fermi-contact 
and dipolar fields B F = 41 T and B dip = 14 T that 
are of opposite sign. As they point out, the value for 
B-p should not vary much for the same ion in compara- 
ble coordinations, quoting values for Mn 3+ in Ti02 and 
in MnFc 2 4 of, respectively, B F = 42 T and 36 T (with 
-Sdip = 12 T and 11.5 T, again of opposite sign). We may 
therefore assume similar values for these fields in Mn6. 
As is often done, in order to estimate the resulting nu- 
clear energy level splittings responsible for the measured 
hyperfine specific heat, we have in sections III B and D 
approximated the net time average of the diagonal part 
of the (slightly anisotropic 55 ) hyperfine interaction by an 
effective isotropic scalar interaction Hm ,i s = ^4s • I, where 
A = — 7 n ftBhf/s, s = 2 and A/k-Q = 8.7 mK (correspond- 
ing to i?hf = 34 T). As for the other nuclei present in 
Mn6, the only sizable contribution to be expected is that 
arising from the long-range dipolar interactions of the 
proton spins with the surrounding electronic spins. Since 
these dipolar fields are small, their contribution may be 
neglected in zero applied field as compared to that of the 
55 Mn nuclei, and only becomes substantial for applied 
fields exceeding 1 T. Although transferred hyperfine in- 
teractions with nuclei on the oxygen and Br ions that 
are directly bonded covalently to the Mn atoms could be 
substantial^ we may nevertheless neglect their contri- 
butions in comparison with the other ones in view of the 
low abundance of the 17 O isotope and the low number of 
Br atoms present. 

Considering next the nuclear spin relaxation (NSR) we 
remark that the longitudinal NSR rate, 1/T± n , is given 
quite generally by the expression 



1 



E„ 



7E3 



2 

m ' 



(10) 



where W m , n denotes the probability for a transition be- 
tween nuclear energy levels m, n induced by the pertur- 
bation considered. As mentioned, we assume the main 



source for NSR to be the time-dependent fluctuations of 
the electron spins s(t) that produce fluctuating compo- 
nents h(t) of the hyperfine field. The theory has been de- 
veloped by Moriyapl2i£ii§£ on basis of the general theory 
of magnetic resonance absorption of Kubo and TomitaiS 3 . 
Two possible sources for the electron spin fluctuations 
have to be considered, namely electron spin-lattice re- 
laxation, characterized by the longitudinal electron spin- 
lattice relaxation rate 1/Ti e , and spin-spin relaxation, 
due to magnetic interactions (dipolar or exchange) be- 
tween the electron spins and characterized by the trans- 
verse electron spin relaxation rate 1/T2 C . Obviously, hy- 
perfine interaction terms producing NSR should involve 
operator combinations as I ± s z , I + s~ and I~s + , I + s + 
and I~s~. The first of these distinguishes itself in that a 
nuclear spin flip is not combined with an electronic spin 
flip. It follows that this first type involves transitions at 
NMR frequencies u> = uj n , the others at ESR frequencies 
lu = uj c ±uj n w lu c . In the high-T approximation, the NSR 
rate can be expressed in terms of the spectral densities 
/j Q of the two-spin (i j) and autocorrelation (i = j) 
functions (sf (0)s"(i)) at these frequencies as 

^ = |^ + i)E[ A ^//(^) + B ^fM' ( n ) 



/j» = / (s?(0)s«(t))e— t dt(a = ±,z). (12) 

J — oc 

The coefficients Aj and Bj are constants depending on 
the details of the hyperfine interactions. For the ease of 
discussion we approximate the Fermi-contact interaction 
by an isotropic scalar on-site hyperfine coupling and ne- 
glect possible orbital contributions. Adding the dipolar 
interaction, the hyperfine Hamiltonian becomes 



Hm = Ali ■ Si 



^-i ' S j _ Q (Ij ' r ij){ S j ' r ij) 



, (13) 



where j c is the electronic gyromagnetic ratio. For the 
terms responsible for NSR one then obtains 



n'= i?Y, D ^) 



(l/2)Asr +J2(D 



s - 



t) 



+ c.c. 



(14) 



where and denote the components of the dipo- 
lar coupling tensor connecting if with Sj and with s. , 
respectively, and c.c. stands for complex conjugates. It 
is important to note that for a pure scalar hyperfine in- 
teraction only the transverse spectral densities f^(uj e ) 
do appear. Assuming an exponential decay of the spin- 
correlation functions, the scalar interaction leads to the 
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1 



— = -A 2 / (s+(Q)s-(t))e^-^ t dt 
3 V ' 1 + (w e - w n ) 2 T 2 



(15) 



As noted above, relaxation then requires energies fiw e 
of the order of the electronic level splittings. By con- 
trast the dipolar interaction contains terms of different 
symmetryjSiSi so that it contributes to both transverse 
and longitudinal terms. In particular it contains the op- 
erator Df- oc — (3/2) sin#cos#e~^ that may induce a 
nuclear flip unaccompanied by an electron flip, thus in- 
volving the much smaller energy Huj n . The NSR rate due 
to this process, after averaging over the angular depen- 
dence, is obtained asS 



1 _ 3 /Mo 



WW) 



Considering now first the case that the electron spin 
fluctuations arise from electron spin-lattice relaxation, 
and assuming again an exponential decay of the autocor- 
relation function, in this case with longitudinal relaxation 
rate 1/T le : 

(s z (0)s z (t)} = i s ( s + l)exp(-VT le ), (17) 
the NSR rate by this process is found to be given by 
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(18) 



an expression first derived by Bloembergen^i for the nu- 
clear spin-lattice relaxation by paramagnetic impurities 
in diamagnetic crystals. The superscript (EZ) is added 
to indicate that this nuclear relaxation is driven by the 
spin-lattice relaxation of the Electron-Zeeman reservoir. 
Comparing Eqs. l|15fl and i|18fl it is clear that, unless the 
electronic linewidths would be comparable to the level 
splittings lu c , the latter process will outweigh the previ- 
ous one by the large factor (w c /aj n ) 2 . In most cases of 
interest one further has w n Ti c 3> 1, so that one may write 
approximately 
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(19) 



Here £?di P = {^o/^)^lc\/s{s+T)r~ 3 stands for the 
electronic dipolar field at the nuclear site, and i?tot = 
B a + Bm = w n /7 n is the total field responsible for the 
nuclear Zeeman splittings^ 

In the next step we have to compare this result with the 
NSR rate arising from spin-spin interactions, which for 



Mng amount to the dipolar interactions between the elec- 
tronic cluster spins S = 12. We note that since the total 
spin S = 12 of the Mng cluster results from the strong fer- 
romagnetic intramolecular exchange between the atomic 
spins s = 2, the fluctuations of the total spin are obvi- 
ously related to those of the constituting atomic spins 
and vice versa. Thus, although both the hyperfine inter- 
actions and the electron spin lattice coupling basically 
involve the atomic spins, the atomic spin fluctuations 
nevertheless are in a one-to-one relationship with those of 
the cluster spins. Due to such spin-spin interactions the 
spectral density will no longer be given by a Lorentzian. 
Instead of the exponential decay of the correlation func- 
tions, Eq. (|T7|) . one usually assumes a Gaussian approx- 
imation for the autocorrelation functions£2i§ii§2: 



= ^s(s + l)exp(-w 1 2 nt i 2 ) 



(16) For the longitudinal NSR rate one obtains: 
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(20) 



(21) 



where (ED) indicates that this process is driven by fluc- 
tuations in the Electron-Dipolar reservoir. For the trans- 
verse relaxation one finds similarly: 



1 



1 



9T (ED) 



[l + exp(-w 2 /2w 2 nt )], 



(22) 



from which it follows that l/Tzn ~ 1/T^ D \ Here cj; nt 
stands for the electronic dipolar spin-spin interaction^ 5 , 
which in our case can be estimated from the dipolar or- 
dering temperature, fiwi n t ~ feeTc an d also corresponds 
to the electronic TSSR rate, ui- m t ~ 1/T 2e . 

At this point it is important to emphasize that the 
above derivations are essentially only valid at high tem- 
peratures and low applied fields, since the effects of po- 
larization of the electronic spins by the applied field have 
been neglected. As noted already by MoriyaS and in 
later work on dynamic polarization^Si the more the 
electron spins become polarized, the less they will be able 
to relax the nuclear spins. To account for this, one should 
replace the time dependencies of the electronic spin s(t) 
by its fluctuating part, 6s(t) = s(t) — So, where Sq denotes 
the thermal average of s. Thus, instead of an expression 
as in Eq. (|17|) for the decay of the electronic spin, one 
should takeSi 

(***(0)to«(t)) = ((*"(o)-*5)(**W-*5)) 

= ((s z - s z ) 2 ) ex P H/T lc ) 

= S(d/dX)4(X)exp(~t/T lc ) 

= S S (d/dX)B s (X)eM-t/Ti c ), (23) 

using the fact that the thermal average of each Mn 3+ spin 
s z is given by sBs(X), where X = gfiBB a S/ftBT, s = 2, 
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and Bs is the Brillouin function for the total molecular 
spin S. Accordingly, the expression ltT7|l for the electron 
spin autocorrelation function should be multiplied by the 
factor 3S(s + 1) dBs/dX. Restricting in what follows 
to the simplest case of spin S = 1/2, as appropriate for 
the present experiments in the high-field / low-T range 
where only the two lowest lying electron Zeeman states 
are relevant, this factor reduces to (1 — tanh 2 X)/2, with 
X = g f i B B a /2k B T. We thus obtain for 1/T^ Z) instead 
of Eq. l(T§|l the relation 
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dip 
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(l-tanh 2 X) (24) 
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whereas instead of Eq. (|21|l one has now: 
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■(1 - tanh 2 X) exp(-w 2 /2u; 2 nt ). (25) 



In both cases, since tanhX gives the degree of polariza- 
tion of the electron spin, one observes that when this ap- 
proaches unity the nuclear relaxation rate goes to zero, 
as to be expected. For the electron-dipolar relaxation 
channel one should notice that, although the actual elec- 
tronic linewidth 1/T% e strongly depends on the electronic 
polarization (the second moment of the absorption line 
is proportional to 1 — tanh 2 X, cf. Ref. I52T) . u>- m t in i|25|) 
is still given by the dipolar coupling as calculated in the 
high-T limits 

Proceeding next to compare the above predictions with 
the high-field NMR experiment, we may already notice 
that Eq. I|25|) yields the right order of magnitude for 
1/T2n- From the value of T c « 0.16 K, we deduce the 
electronic dipolar broadening to be w m t «2x 10 10 rad/s. 
With NMR frequencies of order uj n sa 1.5 x 10 9 rad/s the 
factor exp(— w 2 /2w 1 2 nt ) becomes ~ 1. Further, we have 
Btot « 30 T and B dip a (l/3)B tot » 10 T, yielding 
7n-Bdi P ~ 6 x 10 s rad/s. For applied fields B & > 5 T 
the polarization correction factor (1 — tanh 2 X) becomes 
of order 10~ 3 to 10~ 4 . From Eq. J25J! with the numeri- 
cal factors quoted above we thus find the prediction (cf. 
Fig. HOI dotted line): 1/T 2 ( ™ ) w 10 4 to 10 3 s" 1 for the 
transverse NSR-rate arising from electron spin-spin in- 
teractions, i.e. in the same range as the experimental 
transverse rate. Conversely, the data clearly show that 
1/Tin <C 1/T2 n , contrary to the prediction of Eq. (PSJ. 
Indeed, this process basically only establishes the ther- 
mal equilibrium between the nuclear and the electronic 
spin systems, i.e. without considering the relaxation of 
the latter toward the phonon bath. For the complete de- 
scription of the nuclear spin-lattice relaxation process we 
obviously have to investigate the spin-phonon coupling 
mechanism as well. 

In order to estimate the electronic 1/Ji e , we remark 
that the electron spin-lattice relaxation rate arising from 
transitions between the two lowest Zeeman levels of the 



S = 12 multiplet due to the direct process will be given 
by the sum of the transition rates lUf and due to 
absorption and emission of phonons, respective!^ 

1 . , 

— =w T +w x . (26) 

J-le 

Since the phonon modulation of the crystal field can be 
expected to be the main source of this coupling, we may 
apply to Mn 6 the calculations developed by Leuenberger 
and Loss^ obtaining 

= K- Ph ( ff /W exp(2 ^_ r (27a) 

V e - ph ( 9f iBB a ) 3 - _\ nv , , (27b) 



1 - exp(-2X) 



V f 



c — ph 



D 2 S{2S-l) 2 



(27c) 



where S — 12, D/k B = 0.013 K is the uniaxial anisotropy 
constant, p — 1.45 g/cm 3 is the density and c s the sound 
velocity. Within the Debye model, the latter is obtained 
from the experimental Debye temperature 0d = 29 K as 



fc B 9 D ( 6ttN a \ 



-1/3 



1.3 x 10 3 m/s. 



Substituting into Eq. (|2*F)|) yields 

1/Ti c w 104 S 3 cothX, 



(28) 



(29) 



with B a in Tesla. (It should be noted that the value 
calculated for V c - P h is very sensitive to the values used 
for 0d and D/ks so that it obviously is subject to a large 
uncertainty margin). For instance, B a = 5 T and T = 1 



K yields 1/Ti e 



10 4 



Because of the very small 



value of the anisotropy constant D in Mn@, 1/T le is thus 
expected to be much lower than the typical values ~ 10 7 
s _1 found, e.g., in Mni2-ac. This also implies that a 
model for the TSSR rate based on the random changes in 
local hyperfine field due to electron-phonon excitations, 
as recently used to describe 1/X2n in Mn^-acS would 
lead in this case to a quantitative estimate that is about 
three orders of magnitude lower than our experimental 
result. 

Nuclear relaxation to the lattice can now occur in two 
ways, either directly via the spin-lattice relaxation fluc- 
tuations of the individual electron spins, or in a two-step 
process by spin-spin relaxation to the electron dipolar 
reservoir followed by relaxation to the lattice. The direct 
spin- lattice relaxation (single-ion) process is described by 
Eq. H24f> . which becomes 



1 
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2i I -^dip \ /?•'. ... ,i. v, i i.i 



B u 



B 3 cothX(l-tanh^X). (30) 



For large X, cothX w 1 and (1 - tanh 2 X) w 
4exp(-2X). With B dip /B tot ps 1/3 one obtains 

» 9 Bl cxp(-gp B Bjk B T), (31) 
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1/T 1n « ED > S 




FIG. 11: (Color online) Block diagram of the nuclear and 
electronic spin systems involved in the relaxation process at 
high fields, and the relative rates of energy transfer. 



(with B a in Tesla). As we have seen, both the specific 
heat and the NMR data yield a held dependence of 1 /Ti n 
that is the same as in Eq. QH1JI. with a prefactor of about 
40. The prediction of Eq. I|31[l is therefore qualitatively 
satisfactory but quantitatively slightly too low. 

Next we consider the two-step relaxation process on 
basis of the spin-spin interaction process. Intuitively this 
is easily understood as follows: relaxation of the nuclei by 
spin- spin interactions involves 1/T2 e which will be of or- 
der 10 9 Hz or higher, implying that the electron spin-spin 
interactions can be very effective in relaxing the nuclear 
spins. However, relaxation is then toward the electron 
spin system, and the ultimate relaxation to the lattice 
has to occur in a second step. This situation has often 
been met for nuclear relaxation in magnetic crystals or 
in diamagnetic insulators with paramagnetic impurities, 
notably in connection with the phenomenon of dynamic 
nuclear polarization&S2i2& In the theoretical treatments 
it has been proven necessary to consider the Zeeman term 
and the spin-spin interaction term in the Hamiltonian of 
a spin system (electronic or nuclear) as separate energy 
reservoirs, to each of which separate temperatures can 
be assigned that may diff er q uite substantially from one 
another (see Refs. Ir35l and l70j) . 

Applied to our present problem, this leads to the block 
diagram sketched in Fig. ^2 (the Nuclear-Dipolar reser- 
voir is omitted here since its energy is so small that 
it plays no role at the relevant temperatures). The 
Electron- Zeeman (EZ) and Electron-Dipolar (ED) en- 
ergy reservoirs will be at the same (lattice) tempera- 
ture in zero applied field. However, when with increasing 
field the electronic level splitting ui c starts to exceed the 



electronic dipolar broadening, the two reservoirs become 
progressively separate entities, characterized by different 
temperatures and largely different heat capacities. This 
arises since the EZ reservoir is strongly coupled to the 
phonon bath and can be considered to remain in equilib- 
rium with the lattice regardless of any nuclear relaxation 
event. The ED reservoir on the other hand, has a heat 
capacity that decreases rapidly with held. Although the 
ED reservoir is coupled to the lattice at a rate « 2/Ti c , 
i.e., twice the electron spin-lattice relaxation rate^S it is 
also strongly coupled to the nuclear Zeeman (NZ) system 
at a rate given by Eq. i.e., of the order of 1/T 2n . 

In high fields, therefore, the nuclear relaxation will take 
place in two stages. In the first stage the NZ and ED 
systems will become rapidly in equilibrium at basically 
the same fast rate that determines the experimentally 
observed 1/T2 n . Subsequently, the coupled NZ + ED 
systems will relax toward the lattice at the much slower 
rate 
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(32) 



which for the case Cnz S> Ced would become roughly 
equal to 2/T ic (Ced /Cnz)- Here the symbols Ced and 
Cnz stand for the (field-dependent) specific heats of 
the ED and the NZ reservoirs. The situation is seen 
to be analogous to the phonon-bottleneck phenomenon, 
well known in paramagnetic relaxation. Applied to our 
present case, we calculate Cnz from Eq. [3] to be of order 
0.006 R at T = 0.9 K and B = 5 T (cf. Figs. 6 and 7), 
which value depends only weakly on applied field. Ced 
in zero field can be estimated astS 
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(33) 



after averaging over the angular dependence of the dipo- 
lar coupling, as appropriate for an unoriented pow- 
der sample. Using the same lattice parameters as for 
the Monte Carlo simulations yields Ced ~ 0.004 R 
at T = 0.9 K. Ced then depends on the electronic 
polarization as Ced(A) = Ced(0)(1 — tanh 2 X) rj 
4 Ced(0) exp(— g^ B B a /k B T), since for high fields only 
the lowest Zeeman level is available for the electron spins. 
All this leads to a global rate for the nuclear spin-lattice 
relaxation 



1 



208 Bl coth(A) 



0.004(1 — tanh 2 X) 



0.006 + 0.004(1 



550 Blcxp(-giJ B B/k B T), 



tanh 2 X) 

(34) 



that has the same field dependence as Eq. (|31(l . but 
with a prefactor of about 550 instead of 9, meaning that 
the two-step spin-spin relaxation should be the fastest 
process by almost two orders of magnitude in the high- 
field region. The solid line in Fig. ^| is obtained from 
Eq. (|34|l but assuming a prefactor of order 25. 

At this point it is important to recall that the prefac- 
tors in both Eqs. (|31|l and (|34|l are affected by a large 
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numerical uncertainty originating from the expression for 
1/T le , which contains the elect ron-phonon coupling con- 
stant T4_ph, Eq. H27c|l . This constant is proportional to 
D 2 and O^ 5 , and both these quantities have rather large 
error bars. However, the values of D and 0d influence 
in the same way the relaxation rates via the Electron- 
Zeeman and the Electron-Dipolar channels, thus the lat- 
ter is expected to dominate in any case by almost two 
orders of magnitude. Due to their strong influence on 
1/T lc , allowing both D and 0d to vary by only a factor 
1.5 would already yield the correct quantitative prefactor 
in Eq. J33J|. 

Summarizing the results of this section, we may state 
that both the longitudinal and transverse nuclear relax- 
ation that we observe at high applied fields are in excel- 
lent qualitative and even quantitative agreement with the 
model based on fast dipolar relaxation of the hypcrfinc- 
coupled Nuclear- Zeeman system to the Electron-Dipolar 
reservoir, followed by much slower relaxation of the com- 
bined systems via the electron spin-lattice channel. The 
direct nuclear spin-lattice relaxation process by single- 
ion electron spin-lattice relaxation predicts a similar field 
dependence but is calculated to be much slower at high 
fields. At low fields one will have uj c « a->i nt , so that 
the Electron-Dipolar and the Electron- Zeeman systems 
will become "on speaking terms" , and a subdivision of 
the two electron spin reservoirs is no longer valid. In this 
range, however, NSR by means of the scalar hyperfine in- 
teraction, Eq. I|15[l . should also become important (since 
then no longer uj n <C uj Ci whereas Ti c Ti e ). 



V. CONCLUDING REMARKS 

In conclusion, our experiments on Mng show that 
dipole-dipole interactions between molecular magnetic 
clusters may indeed induce long-range magnetic order at 
low temperatures if the anisotropy is sufficiently small. 
Spin-lattice relaxation is then fast enough to produce 
equilibrium conditions down to the low temperatures 
needed. We should add that similar conditions could in 
principle also be reached in the highly anisotropic clus- 
ter systems, for which it was shown that by applying 
magnetic fields perpendicular to the anisotropy axis, the 
spin-lattice relaxation can be tuned and made similarly 
fast through the process of magnetic quantum tunneling. 
However, it is very rare to observe magnetic ordering phe- 
nomena in those systems. One exception known to date 28 
originates from an unusually high tunneling rate already 
in zero field. Instead, for the majority of the anisotropic 
clusters it is likely that, given the magnitude of the fields 
needed to have a considerable increase of the relaxation 
rate (B± > 1 T), any longitudinal component of the field 
would create a Zeeman splitting that is much larger than 
the energy involved in the magnetic dipolar ordering. We 
found indeed that in Mng the ordering transition is re- 
moved already for relatively small fields (~ 0.5 T). In a 
recent neutron diffraction experiment on a Mni2-ac sin- 



gle crystal, however, Luis et al~ achieved an extremely 
accurate alignment of the field (to within 0.1 degree) and 
obtained evidence for a ferromagnetic phase induced by 
the transverse field. 

We have also studied the nuclear spin dynamics of 
Mn6, both directly by NMR experiments and through 
the hyperfine contribution to the field-dependent spe- 
cific heat. The agreement between the two techniques is 
very good, and also provides an interesting comparison 
with the nuclear spin dynamics in the anisotropic single- 
molecule magnet Mni2-ac. Both qualitatively and quan- 
titatively, the nuclear magnetic relaxation data turn out 
to be in good agreement with predictions obtained from 
theories developed earlier for relaxation in paramagnetic 
crystals and for dynamic polarization. In high fields, the 
observed nuclear relaxation is dominated by electron spin 
fluctuations arising from the dipolar interactions between 
cluster spins. In spite of the large Zeeman splittings be- 
tween the cluster spin levels produced in such high fields, 
these fluctuations are able to relax the nuclei through 
the dipolar part of the hyperfine interaction. In this field 
range the electron dipolar and the electron Zeeman sys- 
tem are basically decoupled. Relaxation of the nuclear 
spins then proceeds in two steps, namely an initial rapid 
relaxation to the electron dipolar system via the electron 
spin-spin interaction channel, followed by a much slower 
relaxation of the combined nuclear-electron spin systems 
to the lattice through the electron spin-lattice channel. It 
is of interest in this regard to note that the values for the 
longitudinal and transverse relaxation rates as observed 
for Mni2-ac in the low-temperature (T < 0.9 K) quan- 
tum regime^ where the electron spin fluctuations can be 
attributed to quantum tunnelling of the cluster spins, fall 
only slightly below the present observations for Mn6 at 
T = 0.9 K. Below the blocking temperature (~ 3 K) the 
cluster spins of Mni2-ac become almost fully polarized 
even in zero field due to the strong crystal field splittings 
of the electron spin levels, the distance between the first 
excited state from the ground state amounting to more 
than 10 K. The temperature independent value found 
for the transverse nuclear relaxation rate, l/T^n ~ 100 
s -1 , could be well explained in terms of intercluster nu- 
clear spin diffusion, i.e., nuclear flip-flops arising from the 
dipolar interaction between nuclear spins in neighboring 
clusters. The same physical mechanism should also put a 
lower bound to the transverse nuclear rate in Mng, which 
is, however, not relevant due to the presence of the faster 
spin-spin relaxation process. 

For the longitudinal rate 1/Ti„ for Mni2-ac a value of 
ss 0.03 s _1 is found below 1 K, slightly depending on 
temperature. This value agrees with the nuclear-spin- 
mediated tunneling rate estimated for the fast-relaxing 
molecular spins in Mni2-ac. As argued by Morello et 
al.r^ the tunneling process can at the same time provide 
a relaxation channel for the nuclei to the electron-dipolar 
system. Similar to the above-discussed case of Mng in 
high field, relaxation to the lattice should then occur in 
a second step through the electron spin-lattice coupling. 
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the grain. In case of a ferromagnetic material, one usu- 
ally only takes the shape-dependent demagnetizing cor- 
rection into account. An argument for this may be found 
in that the magnetization process for the ferromagnet is 
mostly determined by the mobility of the domain walls, 
which will react to the macroscopically averaged internal 
field. For simplicity, we first consider the case of zero 
magnetocrystalline anisotropy for which the demagneti- 
zation factor as well as the magnetization and fields can 
be treated as scalars. We thus obtain for the internal 
field Hi inside the grain 
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Hi — Hi oc — N gla i n M 

= H a - N giain M - fM(N cont - 1/3) (A2) 

From the definition: Hi = H a — N c gM, we thus finally 
find 



APPENDIX A: ESTIMATION OF THE 
EFFECTIVE DEMAGNETIZING FACTOR 

In this Appendix, we derive an approximate expres- 
sion for the effective demagnetizing factor N c s appro- 
priate for the cylindrically shaped container filled with 
the grains. This problem is notoriously complex and can 
only be approximately solved. In a first step the relation 
should be found between the applied field H a and the 
local field H\ oc acting on a reference grain in the con- 
tainer. Approximating the grains as point dipoles, the 
difference (i?i oc — ^a) will be due to the contributions 
to the field arising from all the other dipoles inside the 
container. Adopting the well-known Lorentz construc- 
tion, the dipole summation is split into one inside a (suf- 
ficiently large) sphere around the reference grain and a 
contribution from the dipoles outside this sphere. For 
this second contribution the dipoles are usually assumed 
to form a homogeneous continuum so that it is just pro- 
portional to the difference in demagnetizing factors of 
the container {N coni ) and of the sphere (1/3). As for the 
first summation, it would be zero for a cubic arrangement 
of the dipoles. This will not be the case here since the 
grains are randomly packed, but since a valid estimate is 
not easily obtained, and we may expect it to be small, 
we shall just neglect it. One then obtains 



H a - fM(N cont - 1/3) 



(Al) 



Here / denotes the filling volume fraction of the grains in 
the container and M is the magnetization of the grains. 

In the next step we have to correct H\ oc for the dipolar 
contributions arising from the magnetic material inside 



N cS = A graill + f(N cont 



1/3). 



(A3) 



If, by contrast, the grain has uniaxial magnetic 
anisotropy, we should distinguish between the parallel 
Xi ,i| and the perpendicular Xi,± intrinsic susceptibilities in 
the response to the internal magnetic field. These suscep- 
tibilities depend on the magnitude of the anisotropy pa- 
rameter D. They can be obtained, from the numerically 
calculated eigenstates and eigenvalues of the spin Hamil- 
tonian, u sing the Van Vleck's formalism as described in, 
e.g., Ref. 1x3- I n this case Eq. i|A2[) becomes 



Hi — Hi oc — N„ 



,XiHi, 



(A4) 



where N gla i n and Xi are respectively the diagonal demag- 
netizing and intrinsic susceptibility tensors. By combin- 
ing Eas. HAl|) and (|A4() it is possible to find a relationship 
between the measured susceptibility x an d the two com- 
ponents of Xi- For the case when the anisotropy axes are 
randomly oriented in the sample, we find 



X 



Xcff 



1 + fXes(N cont - 1/3) 



(A5) 



where the susceptibility XcS corrected for the demagne- 
tizing factor of the grains equals 



Xcff 



2X1,. 



Xi 



3(l+Xi,x/3) 3(l+xi,||/3) 



(A6) 



This relationship was used to calculate the theoretical 
powder susceptibilities shown in Fig. [5J 
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